Unidirectional Diagonal Order and 3D Stacking of Charge Stripes in Orthorhombic 

Pr 167 Sro.33Ni0 4 and Ndi. 67 Sr .33NiO 4 . 
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The interplay between crystal symmetry and charge stripe order in Pri.67Sro.33Ni04 and 
Ndi.67Sro.33Ni04 has been studied by means of single crystal x-ray diffraction. In contrast to 
tetragonal Lai.67Sro.33Ni04, these crystals are orthorhombic. The corresponding distortion of the 
MO2 planes is found to dictate the direction of the charge stripes, similar to the case of diagonal 
spin stripes in the insulating phase of La2-iSr :E Cu04. In particular, diagonal stripes seem to always 
run along the short a-axis, which is the direction of the octahedral tilt axis. In contrast, no influ- 
ence of the crystal symmetry on the charge stripe ordering temperature itself was observed, with 
Too ~240 K for La, Pr, and Nd. The coupling between lattice and stripe degrees of freedom allows 
one to produce macroscopic samples with unidirectional stripe order. In samples with stoichiometric 
oxygen content and a hole concentration of exactly 1/3, charge stripes exhibit a staggered stacking 
order with a period of three Ni02 layers, previously only observed with electron microscopy in do- 
mains of mesoscopic dimensions. Remarkably, this stacking order starts to melt about 40 K below 
Too- The melting process can be described by mixing the ground state, which has a 3-layer stacking 
period, with an increasing volume fraction with a 2-layer stacking period. 

PACS numbers: 74.72.Dn, 74.25. Ha, 61.12.-q 
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I. INTRODUCTION 

In recent years the nickelate I^-xSr^NiO^+i has been 
studied intensively because of its similarity to the high 
temperature superconductor I^-xSr^CuO^ 1 ' 2 ' 3 ' 4 ' 5 ' 6 Of 
particular interest is that, in both materials, incommen- 
surable spin correlations and lattice modulations are ob- 
served which are consistent with the concept of charge 
and spin stripes^ Undoped, both systems are insulat- 
ing antiferromagnets (AF). Doping with Sr (x) or O (S) 
introduces hole-like charge carriers, at a concentration 
p = x + 25, into the NIO2 and Cu02 planes, which 
leads to a suppression of the commensurate AF order. 
However, with increasing hole concentration both sys- 
tems show a strong tendency towards a frustrated elec- 
tronic phase separation^ For certain compositions the 
holes segregate into one-dimensional charge stripes, sep- 
arating intermediate spin stripes with low hole concentra- 
tioni In La2_3;Srj;Ni04, stripe order results in an insu- 
lating ground stated Stripes run diagonally to the MO2 
square lattice and are most stable for x = 0.33, with 
the charges ordering at Tco ^240 K and the spins at 
Tqo ^190 Ki2*LiLi£ Stripes can be identified by various 
tcchniquesi^iiiiiSii 3 ^ 1 ^ 5 ^ 1 ^ 1 ^ Here we focus on the char- 
acterization of the charge stripes with x-ray diffraction by 
probing the lattice modulation associated with the spa- 
tial modulation of the charge densityiii 

We are particularly interested in the response of the 



charge stripe order to changes of the crystal lattice sym- 
metry. For certain cuprates it was shown that lattice 
distortions can pin stripes or influence their orientation 
with respect to the lattice. The most prominent exam- 
ple is observed in La2- x Ba x Cu04 and Nd or Eu-dopcd 
La2- a: Sr a ;Cu04 around a hole doping of x = 1/8, where 
stripes are pinned parallel to the Cu02 square lattice as 
a result of a structural transition from orthorhombic to 



tetragonal symmetry. 
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In samples with a fully de- 



veloped static stripe order, superconductivity is strongly 
supprcsscdii&SiLSiiS Another example is observed in or- 
thorhombic La2- a; Sr 2; Cu04 below the metal insulator 
(MI) transition (x < 0.06), i.e., in the short-range or- 
dered spin-stripe phase. The ground state is insulat- 
ing, with spin stripes running diagonally to the square 
latticeiSiSi In this case the coupling to the orthorhom- 
bic distortions is evident from the finding that the spin 
stripes always order parallel to the orthorhombic a- 
axisi^ 5 . Note that, below the MI transition, no diffraction 
evidence for charge stripes has been found. Therefore, 
one has to consider the possibility that here the incom- 
mensurate magnetism might emerge from ground states 
different from stripe order, such as helical spin structures 
or disorder-induced spin topologiesiSiiSLS&SiLiffii^i 

To gain insight into the mechanisms of the charge- 
lattice coupling in the cuprates, we have carried out a 
study of the diagonal stripes in the isostructural nickc- 
lates. The particular question which motivated this work 



is: do lattice distortions in the nickelatcs, which are very 
similar to those in the cuprates, have an impact on the 
stripe orientation? In this context, we consider the crys- 
tal structure of L^-zSr^NiC^+a with Ln=La, Pr, and 
Nd. 

As summarized in Fig, flfa). undoped La2NiC>4 (open 
symbols) undergoes two structural transitions: From 
the high-temperature-tetragonal (HTT) phase to the 
low-temperature-orthorhombic (LTO) phase at Tht — 
750 K, and from the LTO phase to the low-temperature 
less-orthorhombic (LTLO) phase at Tlt — 80 K. (See 
Sec. [H] for definition of phases.) With increasing Sr 
content, both transitions decrease in temperature and, 
above a critical Sr content of x c ~ 0.2, the HTT 
phase is stable for all temperatures i22*Si2liSi2& Hence, 
in Lai.67Sro.33Ni04, stripe order takes place in a tetrag- 
onal lattice, which results in an equal population of do- 
mains with stripes running along the orthorhombic [100] 
and [010] directions'*^ Obviously, such stripe twinning 
complicates any anisotropy study of stripes, such as elec- 
tronic transport parallel and perpendicular to the charge 
stripes or magnetic excitations parallel and perpendicu- 
lar to the spin stripes. 

The substitution of trivalcnt La with the smaller, iso- 
valcnt Pr or Nd causes a significant increase of the 
chemical pressure on the NIO2 planes £2£ As a result, 
the HTT/LTO phase boundary shifts to much higher 
temperatures and Sr concentrations so that ,even for 
x = 0.33, the transition takes place above room tem- 
perature (see Fig.nj. Therefore, the Pr and the Nd-based 
systems allow us to investigate the formation of stripes 
under the influence of orthorhombic strain. The effect of 
the Pr and Nd substitution on both the high tempera- 
ture and the low temperature structural transitions has 
been characterized by several groups, although in some 
studies samples with non-stoichiometric oxygen content 
were investigated^S^ffi In Ref. 0, the charge order in 
Nd2-xSra;Ni04 with x ^ 0.33 was studied with the focus 
on the evolution from stripe to checkerboard order. 

In the present paper we focus on the coupling between 
charge stripes and lattice distortions in the LTO phase 
of Pri.67Si'o.33Ni04 and Ndi.67Sro.33Ni04. In particular, 
we find that stripes always order parallel to the short a- 
axis of the LTO unit cell. This will allow one to produce 
macroscopic samples with unidirectional stripe order by 
detwinning the crystals under unidirectional strain. The 
identified stripe orientation is the same as for the short 
range spin stripe order in lightly doped La 2 _ a: Sr :l; Cu04 
with LTO structure, as will be discussed. The criti- 
cal temperature, Tqo, turns out to be the same as for 
Lai.67Si'o.33Ni04, i.e., no significant dependence on the 
orthorhombic strain is observed. Furthermore, we find a 
strong dependence of the stacking order of stripes along 
the c-axis on the oxygen and hole concentration. No- 
tably, crystals containing excess oxygen exhibit strong 
stacking disorder. Crystals with stoichiometric oxygen 
content show a tendency towards a 3-layer stacking pe- 
riod, resulting in a unit cell enlarged by a factor of 1.5 
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FIG. 1: (Color online) (a) Phase diagram of L^-iSr^NiCU 
with Ln=La, Pr, and Nd as obtained from x-ray powder 
diffraction on polycrystalline specimens i—*— For Pr and Nd 
the HTT/LTO phase boundary is shifted to higher x and T 
with respect to La. Data points marked with an "e" are lin- 
ear extrapolations of measurements up to 800 K. Top: Lattice 
parameters a, b, and c as well as orthorhombic strain at room 
temperature. Data for the single crystals at x — 0.33 are in- 
dicated by triangular symbols, (e) Tilt directions of the NiC>6 
octahedra in the different structural phases. 



along the c-dircction. However, only for a hole concen- 
tration of exactly 1/3 does this 3-layer period lead to a 
well defined superstructure. 



II. CRYSTAL STRUCTURE 

Ln2_ x Sr x Ni04+5, with Ln = La, Pr, or Nd, crystal- 
lizes in the K 2 NiF 4 structure and at high temperature 
is expected to assume tetragonal symmetry, space group 
IA/mmm, regardless of the Sr and O contents^Si! 3 . For 
convenience, it is common to index all phases on the ba- 
sis of the orthorhombic \[2a x \[2b x c supercell, in which 
case the space group of the HTT phase is FA/mmm*^ 
As a function of temperature and doping, a total of 
four structural phases are observed, which can be de- 
scribed by different buckling patterns of the NiOg octa- 
hedral network ^^^ The HTT phase is the only phase 
with untilted octahedra, i.e., the MO2 planes are flat [cf. 



Fig. njc)]. In the LTO phase, space group Bmab, octa- 
hedra rotate about the a-axis, which runs diagonally to 
the MO2 square lattice. In the LTT phase, space group 
P^/ncm, the octahedral tilt axis is parallel to the Ni- 
O bonds, which means that it has rotated by a = 45° 
with respect to the LTO phase. Moreover, its direction 
alternates between [110] and [1-10] from plane to plane. 
The LTLO phase, space group Pccn, is an intermediate 
between LTO and LTT with a < 45°. 



III. EXPERIMENT 



Two single crystals, with compositions 

Pri.67Sr .33NiO4 and Ndi. 6 7Sr .33NiO4 and sizes of 
4.5 mm in diameter and several centimeters in length, 
were grown at the Laboratoire dc Physico-Chimic des 
Solides at Orsay by the travelling-solvent floating-zone 
method. In both cases, crystal growth was performed 
in 1 atm of air. The polycrystals used to map the 
phase diagram in Fig. Q] were synthesized by a standard 
solid-state reaction^^i To remove non-stoichiometric 
oxygen interstitials (excess oxygen), the polycrystals 
were annealed under reducing conditions that depended 
on the Sr content. Powder x-ray diffraction patterns 
of the polycrystals for temperatures up to 800 K were 
taken with a standard laboratory diffractometeriSdi 

Single crystal x-ray diffraction experiments were per- 
formed at beamline X22C of the National Synchrotron 
Light Source (NSLS) at Brookhaven and at beamline 
BW5 of the Hamburg Synchrotron Laboratory (HASY- 
LAB). At BW5, the photon energy was set to 100 keVi& 
At this energy, a sample several millimeters thick can be 
studied in transmission geometry, allowing one to probe 
its bulk. In contrast, at Brookhaven 8.1 keV photons, 
which have a penetration depth on the order of 10 ^m, 
were used, so that, here, samples were studied in re- 
flection geometry. The Pri. 6 7Sr .33NiO4 crystal, with 
a polished and twinned [100]/[010] surface, at first has 
been studied as-grown- in-air at X22C and BW5. Subse- 
quently, it was studied again at X22C, after removing the 
excess oxygen (5 = 0) by Ar-annealing at 900°C for 24 h. 
The Ndi.67Sr .33NiO4 crystal was studied at BW5 after 
being subjected to an identical Ar-annealing. For the 
hard x-ray diffraction experiments no polished surface is 
required. 

The specific heat of the Ar-annealed Pri. 6 7Sr .33NiO4 
crystal was measured with a Physical Property Measure- 
ment System from Quantum Design. The a6-plane re- 
sistivity of a bar-shaped, Ar-annealed Ndi.67Sro.33Ni04 
crystal was measured using the four-probe technique. 



IV. RESULTS 



A. Phase diagram 



To put the results for the two single crystals (x = 0.33) 
into context with Lai.67Sro.33Ni04, we refer back to 
Fig. H[a), the Sr-doping phase diagram for Ln=La, Pr, 
and Nd^^i The shift of the HTT/LTO phase boundary 
to higher temperatures Tht and x with decreasing rare 
earth ionic radius is obvious (La 3+ > Pr 3+ > Nd 3+ ). 
In those cases where Tht > 800 K, the transition tem- 
perature was determined by linearly extrapolating the 
orthorhombic strain to higher temperatures. The low- 
temperature structural transition was observed in all or- 
thorhombic samples. However, it changes qualitatively 
as a function of the Sr content. At low x, the transi- 
tion is of the discontinuous LTO/LTLO type whereas, at 
higher x, it becomes continuous. No reliable information 
was obtained for the crossover concentrations, or whether 
these are different for the La, Pr, and Nd-based systems. 
However, with x approaching x c (Tht = 0), the transi- 
tion becomes very broad and the orthorhombicity barely 
decreases below Tlt- For Nd-based samples, similar re- 
sults were reported in Ref. |4{J. Note that while in the 
Pr-based polycrystal with x = 0.35 a weak LTO/LTLO 
is still visible, it is not observed in the Pr-based single 
crystal with x = 0.33 (see next section). In Figs, ^b) 
and (d) we show the Sr-doping dependence of the lattice 
parameters, and in Fig. ^c) the orthorhombic strain at 
room temperature. 

Corresponding data for the single crystals (triangu- 
lar symbols) are in fair agreement with the polycrystal 
data (cf. Fig. 0). The largest deviations concern the 
absolute values of the lattice parameters a and b. It 
seems that the single crystal data are systematically too 
low by about 0.5%. Differences of this magnitude are 
within the error of our single crystal diffraction experi- 
ment, where a, &, and c were determined from basically 
one reflection each. The temperatures of the HTT/LTO 
transition arc 370 K (315 K) for the surface (bulk) in 
Pri. 6 7Sr .33NiO4 (details in Sec. lIVCT> . and an estimated 
1100 K in Ndi.67Sro.33Ni0 4 . In addition, the Nd-based 
crystal undergoes the LTO/LTLO transition at ~125 K. 

The approximate melting temperatures in Fig. ^a) 
are based on measurements of the surface temperature 
of the melt using an optical pyrometer. These results 
were obtained in earlier crystal growth experiments on 
Ln2_ 2; Sr : i : Ni04+5 with x = 0, 0.25, and 1.5 in which the 
skull melting method was applied. We mention that our 
results in Fig. ^ arc consistent with available data by 
other groups, in particular when taking the effect of oxy- 
gen doping into account ^^^JQ^ However, Fig.^may 
be the most coherent account of Ln2_ 2 ,Sr 2; Ni04+5 with 
stoichiometric oxygen content for x < 0.35 and Ln=La, 
Pr. and Nd. 
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FIG. 3: (Color online) Scans performed on Ndi.67Sro.33NiC>4 
at T — 200 K through possible charge stripe peak positions 
(0, 4 — 2e, 3) and (4 — 2e, 0, 3) of domains A (round symbols) 
and B (quadratic symbols) in the orientation matrix of do- 
main A. Arrows indicate expected peak positions in the case 
that stripes would also run along the 6-axis. Enumerated 
scans in (c,d) are indicated in Fig. [5] The H-signs in (c) and 
(d) indicate the instrumental resolution full width. 



FIG. 2: (Color online) Projection along I of the (h, fe)-zone. 
Domain A: round symbols. Domain B: quadratic symbols. 
Tilt between the two domains exaggerated. Fundamental re- 
flections are black, LTO superlattice peaks grey, LTLO su- 
perlattice peaks green (with dot) and white. The charge 
stripe peaks (CO) are red and blue (with cross), and ap- 
pear at incommensurate positions, displaced by 2e along k 
from Bragg peaks. Reflection conditions for I are indicated 
in the legend. Al depends on the stacking order, which we 
discuss in Sec. II V El The three enumerated scans are those 
inFig.Efod). 



B. Structural transitions in single crystals 

In this and the next section, we focus on the tempera- 
ture dependence of the crystal symmetry and its influence 
on charge stripe order. As mentioned earlier, the Pr and 
Nd-based compounds undergo the HTT/LTO transition 
above room temperature. It is well known that below 
this transition, crystals tend to form twin boundaries in 
the afr-plane. In particular, in space group Bmab, up to 
four twin domains can form, resulting in a corresponding 
manifold of fundamental reflections^ In both crystals 
studied, only two out of these four domains are present, 
which we call domains A and B. 

Figure |21 shows a projection of the (Zifc)-zone along the 
/-direction in the case of the LTLO phase with charge 



stripe order, as we find it in the Ndi.g7Sro.33Ni04 crystal 
at low temperatures. The legend shows the /-conditions 
for the different types of reflections. The fundamental re- 
flections of the domains A and B are indicated by black 
symbols. The angle between the two domains amounts 
to ~0.7° at 200 K (exaggerated in the figure). Super- 
structure reflections in the LTO phase, such as (032), are 
indicated by grey symbols. In the LTLO phase, addi- 
tional reflections appear, such as (110) and (302), which 
are indicated by symbols with dot (green) and open sym- 
bols, respectively. Furthermore, the angle between the 
domains is smaller than in the LTO phase. In the LTT 
phase, the reflections of the two domains are merged into 
single peaks. In our Ndi.67Sr .3 3 NiO4 crystal, the LTT 
phase is not observed down to T = 7 K. 



C. Unidirectional Stripes in the LTO phase 

Now let us turn to the weak reflections which appear 
below Tqo ~ 240 K due to the charge stripe order, in- 
dicated in Fig. [3 by symbols with cross for domain A 
(red) and domain B (blue). In both crystals, as well 
as in each of their two domains, these peaks were always 
found along k. This conclusion is based on scans through 
possible charge-peak positions in both domains, similar 
to those presented in Fig. [31^£ Note that all scans where 
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FIG. 4: Temperature dependence of the normalized intensi- 
ties of (a) the LTO superstructure peak, (b) the LTLO super- 
structure peak and (c) the charge stripe order peak for the 
Ar-annealed Ndi.67Srn.33Ni04 crystal. Data collected trans- 
mission geometry (100 keV). 
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FIG. 5: Temperature dependence of the normalized intensi- 
ties of (a,c) the LTO superstructure peak and (b,d) the charge 
stripe order peak in Pri.67Srn.33Ni04. Data on left side in 
transmission geometry (100 keV) on as-prepared-in-air crys- 
tal, and on right side in reflection geometry (8.1 keV) after 
Ar-annealing. 



performed in the orientation matrix of domain A, so that 
the positions of the peaks of domain B appear slightly 
shifted. The results indicate a lattice modulation along 
the &-axis, which means that in the LTO phase stripes are 
parallel to the short a-axis. Obviously, the orthorhom- 
bic strain in the LTO phase dictates the direction of the 
stripes. No change of the stripe pattern was observed in 
the LTLO phase of Ndi.67Sro.33Ni04, which is not sur- 
prising because in the LTLO phase one still has a < b. 

For the Ndi.67Srn..33Ni0 4 crystal, the temperature de- 
pendence of various superstructure peaks is shown in 
Fig. 0] Below room temperature, the (032) peak in- 
tensity first increases slightly, but starts to decrease be- 
low — 125 K. The decrease of the (032) peak coincides 
with a strong increase of the (330) peak, indicating the 
LTO/LTLO transition.^ Above -125 K, the (330) peak 
is very weak and broad but remains visible up to room 
temperature^ The decrease of the (032) peak in the 
LTLO phase is due to a shift of intensity to the (302) peak 
(white symbols in Fig.[2J). Below the charge stripe tran- 
sition at Tco ^ 240 K, the intensity of the (0, 4 - 2e, 3) 
peak first increases steeply, but tends to saturate below 
200 K. There is no obvious change of the intensity of this 
peak due to the LTO/LTLO transition. 

In Fig. we show corresponding data for 
Pri.67Sr .33NiO 4 . As mentioned before, no LTO/LTLO 



transition was observed in this crystal. The figures on 
the left side show 100 keV x-ray data collected in trans- 
mission geometry on the crystal in the as-prepared-in-air 
state (5 > 0). The (032) reflection strongly decreases 
with increasing temperature, indicating a LTO/HTT 
transition temperature of —315 K. Charge stripe order 
sets in at the same temperature as in Ndi.67Srrj.33Ni04, 
although the transition is somewhat broader [Fig. Etb)]. 
The latter effect is due to the non-stoichiometric oxygen 
and its inhomogeneous distribution in the crystal. The 
influence of excess oxygen will be discussed in more 
detail in Sees. ITVDl and ITVEl 

The figures on the right side of Fig. 03 show data 
collected at 8.1 keV in a reflection geometry after Ar- 
annealing (5 ~ 0). The charge stripe transition now 
is much sharper, though the transition temperature re- 
mains the same. On the other hand, the LTO/HTT 
transition is observed at a significantly higher tempera- 
ture. [Note that in this experiment we have studied the 
equivalent (452) reflection.] In addition to the data in 
Fig. we have collected another set of data (not shown) 
in reflection geometry on the as-prepared crystal. Also 
in this case the LTO/HTT transition occurs at a much 
higher temperature (370 K) than in transmission geom- 
etry. Therefore, we think that Tjjt n ot only depends on 
the hole concentration [cf. Fig.Q{a)], but generally seems 
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FIG. 6: (Color online) Comparison of fc-scans through (h, 4 — 
2e,3) and (h,4 + 2e,3) at 10 K for Pn.erSro.aaNiOi (PSNO) 
(a-c) and Ndi.67Sro.33Ni0 4 (NSNO) (d). Scans in (c) per- 
formed at h — 4 (triangular symbols). All other scans at 
h = 0. Data collected (a) in reflection and (b) in transmission 
on as-prepared sample, (c) in reflection and (d) in transmis- 
sion on Ar-annealed samples. The intensities of the peaks at 
k — 4 + 2e were normalized by same factors as their coun- 
terparts at k = 4 — 2e. The H-shaped sign in (d) indicates 
typical instrumental resolution full width. 
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FIG. 7: (Color online) (a) Comparison of i-scans at 10 K 
through charge peak at (h, 4 — 2e, 3) for Pri.67Sro.33Ni04 and 
Ndi.67Sro.33Ni04. For triangular symbols scans were per- 
formed at ft. = 4. All other scans at h — 0. Typical in- 
strumental resolution full width is indicated by H-sign. (b) 
Incommensurability e as a function of temperature. Error 
for e indicated by vertical bar. Symbols and color code for 
samples and experiments the same as in Fig. El 



to be somewhat higher in the surface layer than in the 
bulk. 



D. In-plane Charge Stripe Distance 

Excess oxygen 5 increases the total hole concentration 
p = x + 2(5 in the NIO2 planes. One way to measure p 
is to measure the in-plane distance between the charge 
stripes via the incommensurability e (see Fig. |2J). The 
property e(p) was studied intensively in recent years and 
it is well known that it is not precisely linear^ However, 
around p = 1/3, one finds that e — p is approximately 
satisfied. Hence, deviations from p = 1/3 can be probed 
quite accurately. A precise way to determine e is to take 
one fourth the distance between the simultaneously mea- 
sured reflections at k = In — 2e and k = 2n + 2e. 

In Fig.0wc show representative scans for as-prcpared- 
in-air and Ar-annealed crystals taken at different pho- 
ton energies, as indicated in the figure. The larger the 
distance between the peaks, the larger p. The dashed 
lines are for p = 1/3. By comparing the 8.1 keV and 
100 keV data in Figs. |^a) and (b), we find that in the 
as-prepared Pri.67Srn.33Ni0 4 crystal the oxygen content 



in the surface is larger than in the bulk. Only after the 
Ar-annealing is the hole concentration in this crystal very 
close to the nominal value p = x = 1/3, suggesting that 
S ~ [Fig. Etc)]. In the Ar-annealed Ndi. 67 Sro.33Ni0 4 
crystal we observe p < 1/3, which we attribute to a Sr 
content slightly below x = 1/3, assuming that 6 ~ 
[Fig- El^d)]. Note that for Ndi.67Srrj.33Ni04 we have mea- 
sured the peak at k = 4 — 2e, only. Another feature 
is that, in the as-prepared Pri.g7Sro.33Ni04 crystal, the 
peaks are asymmetric and obviously broader in k than in 
both Ar-annealed crystals, indicating that excess oxygen 
leads to a enhanced distribution of e. 

In Fig. Ufb) we show the corresponding temperature 
dependencies of e. The symbols are the same as in Fig. El 
At low temperatures one can clearly see the deviations 
from e = 1/3 for as-prepared Pri.67Sro.33Ni04 and Ar- 
annealed Ndi.67Sro.33Ni04. It is believed that the low 
temperature values of e represent the true hole concen- 
tration p. However, close to Tcoi in all measurements 
e gravitates towards a value of 1/3. This lock-in effect 
is well known and indicates that the in-plane stripe dis- 
tance prefers to be three Ni sites (1.5x6), i.e., commen- 
surate with the latticei 1:L i 51 i 52 i 53 i 54 Note that in the case 
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FIG. 8: (Color online) Diagrams indicating the positions of 
the charge stripe peaks (red symbols) as well as correspond- 
ing ordering wave vectors for samples with (a) 3-layer stripe 
stacking period and (b) 2-layer period or strong stacking dis- 
order. Full symbols indicate fundamental reflections, open 
symbols LTO superstructure reflections. In (b) we assume a 
2-layer period similar to Fig. HUT c^ . Additional charge peaks 
would appear for I even, assuming the stacking depicted in 
Fig. 1Mb) ^a 



of the Ar-anncalcd Pri.67Sro.33Ni04 crystal, e ~ 1/3 for 
all temperatures, which indicates that this is in fact the 
only experiment where we were truly looking at a sample 
with p = 1/3. 



E. Three Dimensional Charge Stripe Order 

The most dramatic effect the Ar-annealing had on 
the Pri.67Sro.33NiC>4 crystal was a drastic change of 
the correlations between charge stripes along the c- 
axis. Corresponding /-scans through the charge stripe 
peak at ~ 10 K are shown in Fig. |7fa). In the as- 
prepared crystal, a single narrow peak appears at I = 3, 
which is in accordance with results by other groups on 
La 2 - 3 ;Sr a ,Ni04 + ji :L1 i 55 i 56 After the Ar-annealing, how- 
ever, the peak is split with two pronounced maxima ap- 
pearing at / ~ 2.7 and / ~ 3.3. Since these values are 
close to I = 2 + 2/3 and I = 4 — 2/3, the split most likely 
indicates that a long range stacking order of stripes has 
developed along the c-axis with a period of 3/2 x c, cor- 
responding to a period of three MO2 layers. Figure |S| 
shows the locations of the charge stripe peaks in the 
(Okl)-zone as observed for the Ar-annealed (left) and the 
as-prepared Pri. 6 7Sr .33NiO4 crystal (right). The corre- 
sponding ordering wave vectors are gco = (0, 2e, 2/3) 
and g co = (0' 2e, —2/3) as well as gco = (0, 2e, 1), re- 
spectively. Details will be discussed in Sec. |VJ In the 
Ndi.67Sro.33Ni04 crystal, a similarly clear splitting in / 
was not observed, but, as one can see in Fig. |7{a), the 
scan is extremely broad and shows two shoulders at com- 
parable /-positions, plus a central maximum, indicating 
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FIG. 9: (Color online) (a) /-scans through (4, 4 - 2e, 3) of 
Ar-annealed Pri.67Sro.33Ni04 for different temperatures. In- 
strumental resolution full width is indicated by small H-sign. 
(b) Peak positions as a function of T from fit with three peaks. 
The position of the central peak was kept fixed at I — 3. 



the presence of incommensurate peaks obscured by dis- 
order. 

In Fig. E2a) we show the temperature dependence of 
/-scans through (4,4 - 2e, 3) for Pri. 67 Sr .33NiO 4 . To 
extract the peak positions we have fit the data with three 
lines, keeping the position of the central line fixed at 
/ = 3. As one can see in Fig. |H{b), the splitting in / is 
about constant up to 200 K and then starts to decrease. 
Eventually, for T > 228 K the three peaks have merged 
into a single reflection at I = 3, indicating the loss, or a 
qualitative change, of the long-range correlations along 
the c-axis. 

In Fig. OH peaks appear at slightly asymmetric posi- 
tions with respect to / = 3. We have performed ad- 
ditional scans through charge and nearby Bragg peaks, 
which confirm that the peak positions are symmetric. 
The deviations in Fig. [§] arc due to a small misalignment 
of the crystal, as well as the temperature dependence of 
the lattice parameter c. In our data analysis in Sec. [V] 
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FIG. 10: Stripe stacking patterns, (a) 3-layer period com- 
parable to cubic stacking, (b) 2-layer period comparable to 
hexagonal stacking, (c) 2-layer period due to body centered 
stacking with alternating layers of Ni-O-Ni-bond centered 
(black symbols) and Ni-site centered (grey symbols) stripes. 



this deviation is compensated by a small offset in I. 



V. DISCUSSION 

A. Modelling the stacking of stripes 

Previous studies of L^-zSr^NiC^+s have shown that 
the stacking of stripes along the c-axis is mainly con- 
trolled by two mechanisms^*^: the minimization of the 
Coulomb repulsion between the charge stripes, and the 
interaction of the charge stripes with the underlying lat- 
tice. The Coulomb interaction favors a body centered- 
type stacking to maximize the distance between near- 
est neighbor stripes in adjacent planes. On the other 
hand, the interaction with the lattice favors shifts by in- 
crements of half of the in-plane lattice parameter. This 
means that nearest-neighbor stripes do not stack on top 
of each other. For a hole content of p = 0.25 and 0.5, 
a body-centered stacking also satisfies commensurability 
with the lattice. This is not the case for p = 1/3. Here, 
a body-centered stacking is achieved only in the case 
of an alternation of layers with site-centered and bond- 
centered stripes, similar to the stacking shown Fig. llOf c). 
However, our results for Ar-annealed Pri.67Sro.33Ni04 
strongly suggest a stripe ground state with a 3-layer pe- 
riod, consistent with a stacking similar to Fig. HOf a). 
Stripes successively shift by increments of 6/2, i.e., there 
are three different types of planes with stripes at in-plane 
positions 0, 6/2 and b. A stacking pattern such as that 
in Fig. HOf a - ) but with all stripes centered on Ni-sites is 
possible, too. 

This problem can be projected onto that of the stack- 
ing of close-packed layers of atoms, where layers can be 
shifted by 1/3 of the primitive translation, again result- 
ing in three types of planes: A, B and C^£ The two 
most simple arrangements of these planes are the hexag- 
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FIG. 11: (Color online) (a) Calculated scattering intensity 
for /-scans as a function of the probability z, where z = for 
hexagonal close packing and z — oo for cubic close packing 
(after Ref. I57B . The curve for z — 1 (dashed) represents 
complete disorder. The curve for z = 3.5 represents a system 
with a fairly developed 3-layer stacking period, (b) /-positions 
of the maxima as a function of the volume fraction of ABCA 
stacking segments. 



onal packing ...ABAB...^ and the cubic packing ...AB- 
CABC..., corresponding to stripe stacking patterns with 
2-layer and 3-layer periods, respectively [cf. Fig. HUT a.b)]. 
Energetic differences between these two stacking types 
appear only if second-nearest-neighbor interactions are 
included. 

In Fig. lllf a) we show the /-dependence of the scattered 
x-ray intensity for close-packed structures with second 
neighbor interaction, following the calculations by Hen- 
dricks and Teller M£L Stacking faults can be introduced 
by tuning the probability for second neighbor layers to 
be alike or unlike. For details, see section 6 of Ref. |&2l 
The curves in Fig. lllf a) were calculated for different de- 
grees of stacking disorder z, ranging from a predominant 
2-layer period for z -c 1 (central peaks at integer /) to 
a predominant 3-layer period for z» 1 (split peaks at 
J ±2/3 with / even). The dashed line (red) line is for com- 
plete disorder (z = 1), which is characterized by a cen- 
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FIG. 12: (Color online) /-scans through (4, 4 — 2e, 3) charge 
peak of Ar-annealed Pri.67Sro.33Ni04 plus fits (red solid line) 
as described in the text. The dash-dot line (blue) represents 
the structure factor, the dotted line (green) a T-independent 
constant background. Instrumental resolution full width is 
indicated by H-sign in (a). 



tral peak for odd I and no obvious remains of the peaks 
at even ImH In Fig. Illf b) we show that, as the volume 
fraction with a 3-layer period decreases, the splitting in 
I decreases and eventually disappears close to the point 
where complete disorder is reached. 

The experimental results for Ar-annealed 
Pri.67Sro.33M04, in Fig. clearly show that for 
p = 1/3 a ground state with a 3-layer stacking period 
is favored over a 2-layer period. First, this indicates 
that stripes prefer to be all of the same type, i.e., either 
bond- or site-centered, but not mixed as in Fig. HOf c). 
Second, it shows that the Coulomb repulsion between 
stripes in second neighbor planes is significant. The fact 
that in Fig. 02a) at 10 K the splitting in I is somewhat 
smaller than 2/3, indicates a certain degree of disorder 
[cf. Fig. Illf b)]. Furthermore, it is obvious that the 
peak intensity at approximately (4,4 — 2e,2 + 2/3) is 
significantly lower than at (4, 4 — 2e, 4 — 2/3) . We assume 
that this difference is largely due to the structure factor, 
since the Bragg intensity at (442) is about a factor forty 
smaller than at (444) ■££ From the width of the larger of 
the two peaks we have estimated a correlation length of 
£ c = 17.5 A, which corresponds to approximately three 
Ni0 2 layers^i 

In Fig. E| we present fits to the data obtained by the 
matrix method used by Hendricks and Teller^ multi- 



plied by a slowly varying structure factor F 2 for which 
we assume the phenomenological, Gaussian /-dependence 
indicated in Fig. lT2T a). The only variables are the stack- 
ing parameter z and the amplitude. A finite background 
estimated from fc-scans, a small offset in / of -0.04 (expla- 
nation in Sec. 11V EJl , and the structure factor were kept 
constant for all temperatures. Obviously, the model pro- 
vides a good description of the data. Moreover, the pa- 
rameter z gives us access to the temperature dependence 
of the stacking disorder. Note that fits with a different 
F 2 than the one shown in Fig. ll2f a). e.g., with the max- 
imum centered at I = 4, had virtually no effect on the 
T-dependence of z. 



B. The charge stripe melting process 

In Fig. El wc compare the fitted z(T) with var- 
ious properties of the Ar-annealed Pri.67Sro.33NiC>4 
and Ndi.67Sr .33NiO4 crystals, such as the /-integrated 
charge-peak intensity, the inverse in-plane correlation 
length ^j/ 1 of the charge stripe order, the incommensura- 
bility e, the specific heat c p , and the in-plane resistivity 
p a b- A pronounced maximum is observed in c p at 238 K, 
which we identify with Tco, the onset of static charge 
stripe order. Above this temperature, charge peaks be- 
come very small, and £r starts to diverge, i.e., the in- 
plane stripe order melts^ In contrast, the parameter z 
already starts to decrease at ^200 K and reaches com- 
plete disorder (z = 1) at about 228 K. From this we 
conclude that the melting of the 3D charge stripe lattice 
sets in with a melting of the interlayer correlations, leav- 
ing the 2D correlations intact, until these eventually melt 
at Tco as well. 

In Fig. I13f b) one can see that z goes through a min- 
imum and then increases for T > 235 K. This anew in- 
crease goes along with the increase of t;^ 1 and reflects 
the drastic decrease of the stripe correlations above Tco ■ 
Note that, in Lai.67Si'o.33Ni04, a clear minimum in £ t ~ 1 
was observed at about the same temperature^ We be- 
lieve that there the increase of £~ at temperatures below 
the minimum is a precursor of the formation of the 3-layer 
period stacking order. According to Fig. IT^T b). z < 1 at 
the temperature of the minimum, which, if taken seri- 
ous, implies a slight tendency towards a 2-layer stacking 
period. It is certainly possible that, once the stacking 
order starts to melt, stripes can arrange more freely, and 
may assume a body-centered stacking order, for which 
the ordering wave vector gco is the same as in the case 
of disorder [cf. Figs. E^b) and llOf c')]. However, further 
work is necessary to verify if our analysis, which is based 
on the model by Hendricks and Teller, correctly describes 
even such small details. 

The spin stripe ordering temperature Tgo of our 
crystals has not yet been determined. While in 
Lai.67Sro.33Ni04 both Tco an d Tso can be identified by 
weak anomalies in the static magnetic susceptibility of 
the Ni02 planes, this is impossible in the case of our 
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FIG. 13: (Color online) Top: temperature dependence ol var- 
ious properties ol Ar-annealed Pri.67Sro.33Ni04 close to the 
charge stripe order transition, (a) /-integrated intensity ol 
(4, 4 — 2e, 3) charge peak, (b) parameter z from fits ol /-scans 
(see Fig. 1121 , (c) inverse in-plane correlation length from k- 
scans (d) incommensurability e (see Fig. 0, (e) measured 
specific heat (inset) and specific heat anomaly obtained by 
subtraction of 4th-order polynomial; solid line (red) in inset. 
Bottom: afo-plane resistivity and its second derivative ol the 
Ar-annealed Ndi.67Sro.33Ni04 crystal. 



crystals, because of the large paramagnetic contribution 
of the Pr 3+ and Nd 3+ ionsS However, results from neu- 
tron diffraction in Ref.0show that, in Ndi.67Sr .33NiO4, 
spin stripes order at essentially the same temperature as 
in Lai.67Srg 33MO4, with Tgo — 190 K. This means that 
the spin stripe order disappears just about the tempera- 
ture where the c-axis stacking correlations of the charge 



stripes start to melt. In this context it is worth noting 
that the low temperature side of the anomaly in c p has a 
shoulder which extends down to Tso- Measurements of 
c p in Ref. |62| for a Lai.67Sro.33Ni04 single crystal show 
this even more clearly. We believe that the shoulder in- 
dicates that the melting of the stripe stacking order is 
associated with a significant entropy changed 

Fingerprints of the melting of the intcrlaycr stripe or- 
der well below Tqo are also observed in the electronic 
transport. In Ndi.67Sro.33Ni04, in Fig.^f), the onset of 
the charge stripe order at Tqq is announced by a kink in 
log(/O a b), as is also evident from a peak in d 2 \og(p a b) / dT 2 . 
However, the entire transition in log(p a h) is not sharp, 
but extends down to about 200 K, similar to the behav- 
ior of z, Ac p , and the charge-peak intensity. One might 
think of this as a manifestation of the intinrate connec- 
tion between charge and spin stripe order, since it implies 
that the charge order has to be fairly progressed, before 
spin stripe order can occur. 

Another interesting feature of the anomaly in c p is 
its extension towards high temperatures. In this region 
charge peaks become extremely weak and broad. It is 
reasonable to consider that both effects are related to 
slowly fluctuating charge stripe correlations^ Intensity 
from these correlations is picked-up due to the poor en- 
ergy resolution of the x-ray diffraction experiment. An 
inelastic neutron scattering study of the incommensu- 
rable magnetic fluctuations in Lai. 7 25Sr .275NiO4 gives 
indirect evidence for a liquid phase of charge stripes 
above Tco^ 

The present observations are consistent with a stripe 
smectic or nematic phase at T > Tco<^ (Note that this 
assignment is different from that originally made by Lee 
and CheongiSi) 



C. Comparison with TEM data 

Our results concerning the stacking of stripes are 
in qualitative agreement with a recent transmission- 
electron-microscopy (TEM) study on tetragonal 
Lai.725Srrj.275Ni04 with x=0.275ii^ In this study, it 
was shown that, in domains of mesoscopic dimension, 
charge stripes are indeed one dimensional. However, 
the Sr content of x = 0.275 gave rise to a mixing of 
features expected either for x = 0.25 or x = 0.33, as 
well as features unique to this intermediate doping. In 
particular, the obtained average in-plane stripe distance 
(1.75 x b) results in a mixture of site- and bond-centered 
stripes within the same MO2 plane. As a consequence, 
the stacking of stripes along c is strongly disordered, 
with both simple body-centered arrangements as well as 
staggered shifts by 0.5 x b prevalent^ 

Here, we have shown that the orthorhombic symmetry 
of Pri.67Sr .33NiO4 and Ndi.67Sr .33NiO4 can be used as 
a tool to obtain large single crystals with unidirectional 
stripe order, suited to the study of macroscopic proper- 
ties. For x = 0.33 and 5 = 0, stripes in a single MO2 



11 



plane are likely to be either all bond-centered or all site- 
centered, since the stripe distance of 1.5 x b is commen- 
surate with the lattice. In the case of long range in-plane 
order, the Coulomb interaction between stripes in first 
and second neighbor planes becomes crucial and stabi- 
lizes a stacking such as depicted in Fig. HOf a). Thus, 
for the first time, a stripe order with a 3-layer period 
has been observed in the x-ray diffraction pattern of a 
macroscopic crystal. 

It is worthwhile mentioning that there are several stud- 
ies of Lai.67Sro.33Ni04 with e matching the critical value 
of 1/3 quite we ll,S^iSS^ However, a stripe stacking or- 
der of the kind identified here in Pri. 6 7Sr .33NiO4 and 
Ndi.67Sro.33Ni04 has not been reported. These observa- 
tions suggest that the stabilization of the 3-layer stacking 
period may benefit from the fact that, in orthorhombic 
crystals, the in-plane stripe order is unidirectional. 



D. Comparison with La2- I Sr I Cu04 

Besides the obvious differences, the electronic phase 
diagrams of nickelates and cuprates show some interest- 
ing similarities. In the present context it is particularly 
remarkable that, in La2- 2: Sr. E Cu04, as long as x is be- 
low the MI transition, a static diagonal spin stripe or- 
der forms ^&22i2ii It remains an open experimental ques- 
tion, whether corresponding samples also exhibit diag- 
onal charge stripes. In this respect, it is an interest- 
ing finding that, in nickelates with the LTO structure, 
charge stripes prefer to be parallel to the octahedral tilt 
axis, similar to the spin stripes in the cuprates. These 
findings thus provide motivation for further experimental 
searches for charge stripes in the insulating cuprates. 

There are certainly alternative models for the lightly- 
doped cuprates that do not involve stripes 27 ] 28 i 29 i 30 i 31 ; 
however, there is also additional evidence for the cuprates 
suggesting a common coupling of charge and spin stripes 
to an orthorhombic lattice distortion. In I^CuG^+a and 
Lai.88Si'o.i2Cu04, both systems with parallel stripes and 
LTO structure, it was observed that the direction of the 
spin stripes is not perfectly parallel to the Cu-0 bonds, 
but slightly inclined towards the diagonally running oc- 



tahedral tilt axis (|j a), by an angle too large to be ex- 
plained with the orthorhombic lattice distortioniSSiS 7 . Un- 
fortunately, in these compounds, no c harg e peaks were 
detected. However, in Ref. |63 and Ref. |63 it was shown 
that in orthorhombic Lai.875Bao.i25-i;Sr 2 ;Cu04 with x = 
0.075, both spin stripes and charge stripes show the same 
inclination towards the octahedral tilt axis. This strongly 
suggests that, in both classes of materials, nickelates and 
cuprates, the response of charge stripes to orthorhombic 
lattice distortions is similar. 

VI. CONCLUSION 



We have presented a detailed x-ray diffraction study 
which sheds light on the coupling between the charge 
stripes and structural distortions in Pri.67Sro.33NiC>4 
and Ndi.67Sro.33Ni04 single crystals. In contrast to the 
sister compound Lai.67Sro.33Ni04, both crystals undergo 
a HTT^LTO transition well above room temperature, 
so that stripe order at Tqo takes place in an anisotropic 
environment. We find that the orthorhombic strain 
causes the stripes to align parallel to the short a-axis, 
which is also the direction of the NiC>6 octahedral tilt 
axis. In addition to these in-plane correlations, we have 
observed correlations between the planes, which arc 
consistent with a stacking period of three MO2 layers. 
This stacking order is extremely vulnerable to intersti- 
tial oxygen impurities and deviations of the total hole 
concentration from p = 1/3; as a result, it was observed 
in Ar-annealed samples, only. Further, we find that the 
melting of the static charge stripe order is a two-step 
process, which starts at 200 K with the melting of the 
intcrlayer correlations, and is completed at Tqo with 
the melting of the in-plane correlations. Implications 
for the stripe order in insulating La2_ 2: Sr, c Cu04 have 
been discussed. The observation of unidirectional stripes 
in Pri.67Sro.33NiC>4 and Ndi.67Sro.33Ni04 opens up 
the unique possibility to characterize their anisotropic 
properties. 

The work at Brookhaven was supported by the Office 
of Science, U.S. Department of Energy under Contract 
No. DE-AC02-98CH10886. 
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